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Anomalous magnetothermal effects are discussed in the spin-1/2 Heisenberg chain. The en- 
ergy current is related to one of the non-trivial conserved quantities underlying integrability 
and therefore both the diagonal and off diagonal dynamical correlations of spin and energy 
current diverge. The energy-energy and spin-energy current correlations at finite temperatures 
are exactly calculated by a lattice path integral formulation. The low-temperature behavior of 
the thermomagnetic (magnetic Seebeck) coefficient is also discussed. Due to effects of strong 
correlations, we observe the magnetic Seebeck coefficient changes sign at certain interaction 
strengths and magnetic fields. 
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In the last two decades, ID strongly correlated systems 
have attracted immense theoretical and experimental in- 
terest. One of the reasons stems from their unusual static 
and dynamical properties peculiar to ID systems. The 
spin-1/2 Heisenberg chain is one of the most fundamen- 
tal solvable models for ID magnetic insulators and has 
served as a testing ground for many approaches. 

Recently, transport properties of low-dimensional 
strongly correlated quantum systems have been exten- 
sively studied from both theoretical and experimen- 
tal sides. 1 Among them anomalously enhanced ther- 
mal conductivity 2-4 and unconventional large spin dif- 
fusion constants 5 have been reported in experiments on 
one- or quasi ID materials with weak interchain interac- 
tions. These observations indicate the existence of non- 
dissipative transport properties in ID quantum systems. 
Theoretically, the existence of such anomalous proper- 
ties has also been pointed out especially in quantum 
integrable systems. 6 One of the criteria for anomalous 
transport is the existence of a non-zero Drude weight. 

In particular for the spin-1/2 XXZ chain, the Drude 
weight for the spin transport has been obtained by Zo- 
tos for finite temperature at zero magnetic field. 7 In 
the massless regime, the Drude weight is non-zero and 
hence the spin transport is non-dissipative. In the mas- 
sive regime without magnetic field a similar treatment 
yields zero Drude weight at any temperature implying 
dissipative spin transport. This feature, however, contra- 
dicts with numerical or field theoretical approaches and 
is still debated. 8-10 

On the other hand, the thermal Drude weight D t h(T) 
for zero magnetic field has been recently obtained by the 
Bethe ansatz showing that the thermal transport is non- 
dissipative in both massless and massive regimes. 11, 12 At 
low-temperature D t h(T) behaves as D t h{T) = itvfT/3 
where i>f is the velocity of excitations. This universal 
behavior was also found in general systems in which the 
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low-energy excitations are described by a c = 1 confor- 
mal field theory. 13, 14 In the massive regime, we find that 
Ah C 7 ) ~ exp(-S/T)/VT where 5 is the one-spinon (re- 
spectively one-magnon) excitation gap for the antiferro- 
magnetic (respectively ferromagnetic) regime. 15 

In this article we will mainly discuss the thermal trans- 
port and magnetothermal effect in presence of finite mag- 
netic fields. Because the spin-reversal symmetry van- 
ishes for finite fields, the magnetothermal effect does 
appear. 6 ' 16 ' 17 Consequently the thermomagnetic (mag- 
netic Seebeck) coefficient is finite. Moreover we observe 
the magnetic Seebeck coefficient changes sign for certain 
interaction strengths and magnetic fields, which can be 
interpreted as effects of strong correlations. 

Let us consider the spin-1/2 Heisenberg XXZ chain 
with magnetic fields h: 

h 



ft - h kk+i - o y^°fc> 



k=l 



k=l 



h kk+1 = J <{ a+a k+1 + cr k+1 a k 



Here we restrict ourselves to the critical antiferromag- 
netic regime < A < 1 and J > 0. In this case, the 
anisotropy parameter A is conveniently parametrized by 
A = cos 7 (0 < 7 < 7r). 

The spin and energy current operator of the present 
system are written as 

L L 

J s = J }fot a k+i +h.c), Je = iJ }Jhk-ik,hkk+i], 

k=l k=l 

respectively. The transport coefficients are determined by 
the Kubo formula in terms of the correlation functions 
of the above defined current operators. For instance, the 
thermal conductivity is determined by 18 



Re k(uj) = ^Re j L QQ - 



irD th (T)5(uj) + K Te 
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where 



= / dte-' 1 ^ / dT(Ji(-t-iT)Jj 
Jo Jo 



}, = {Q,s}. 



Here the thermal current J7q should be Jq = Je — hj s . 

As already mentioned by Zotos in ref. 6, the energy 
current Je is a constant of motion {[H, Je\ = 0) and 
therefore the thermal Drude weight is finite (D t h(r) > 0) 
at all finite temperatures. Hence the thermal transport 
is non-dissipative. 

For zero magnetic field h = 0, the magnetothermal 
effect is always zero because the system exhibits the spin- 
reversal symmetry. In contrast, for finite magnetic fields 
h > we see {JeJs) > 0, which leads to the diverging 
off diagonal dynamical correlations: 

ReL Qs = 7r(P{J E J s ) - hD s (T))5(u), L sQ = L Qs . 

In this case the thermal Drude weight 

; (JEJs) 2 



D th (T) = 2 (j£)-f3 



D S (T) 



and the magnetic Seebeck coefficient 



S(T) = Rc 



1 Lq s \ _ 1 f (JeJ s ) 1 



TirD s (T)j T\D S (T)T 



(1) 



(2) 



are both finite at finite temperatures. 

To obtain them, we have to evaluate the correlations 
of Je and J s . Let us consider the autocorrelation {J^) 
first. To evaluate this quantity, we introduce the follow- 
ing extended Hamiltonian including as a perturba- 
tion; TL := XqHq — hM + XiJe, where Ho is the Hamil- 
tonian without the Zeeman term hM. The parameters Ao 
and Ai are introduced for later convenience and should 
be taken to be A = 1 and Ai = after taking all 
necessary derivatives. Introducing the partition function 
Z = Tre~P n , we easily see that (J^) is evaluated by 
taking the second logarithmic derivative with respect to 
Ai: (Ji) = d 2 Xi In Z\ Xi=Q /(Lf3 2 ). Note that from sym- 
metric arguments we can prove (J7e) = 0. Applying a 
lattice path integral formulation, we introduce the quan- 
tum transfer matrix (QTM) in the imaginary time di- 
rection. In this formalism Z in the thermodynamic limit 
L — > oo can be expressed as the largest eigenvalue of the 
QTM A(A ,Ai), namely lim L _ 00 (lnZ)/L = lnA(A ,Ai): 

lnA(A ,Ai) = ^+ I a 1 (x+i)ln(l + rr 1 (> + i))dx, (3) 
^ Jc 

where the contour C encloses the real axis (for instance 
we take Imi = i (Imi = — i) for the upper (lower) con- 
tour) and a n (x) = 7sinn7/(27r(cosh7a; — cos 717)). The 
unknown function rj(x) is determined by the following 
non- linear integral equation (NLIE): 19 

In r](x) = p(\ + XtAd^x) + 0h 

+ j> a 2 (x-y-i)ln(l+? ? - 1 (y + i))dy, (4) 

where e(x) = — 2 / nAa 1 {x) and A = 2Jsin7/7. Conse- 
quently we obtain 

(j|) = ^l i lnA(l,A 1 )| Ai=0 . (5) 



In Fig. 1, the temperature dependence of the correlation 
(Je) f° r the isotropic point A = 1 is depicted for various 
magnetic fields. For < < h c (h c = 4 is the critical 
magnetic field for A = 1), wc find that (J^) is linear 
in T at low-temperature T « 1, On the other hand, 
{J%) exponentially decays for h > h c due to the mass 
gap. At high temperature, converges to a constant 
J 4 (l + 2A)/2 not depending on the magnetic field h. 
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Fig. 1. Illustration of {J^) for the isotropic case A = 1 as a 
function of temperature for various magnetic fields h. 



Next we consider the correlation (JeJs)- Since the spin 
current J s is not a constant of motion (except for the XY 
model A = 0), the quantity JeJ% is no longer conserved, 
namely [H, Je<J s ] 0- Hence the above procedure, which 
is useful to calculate the conserved quantities, is not ap- 
plicable directly. Alternatively we use the following non- 
trivial identity relating (JeJs) to {J%), which is valid 
for the thermodynamic limit L — > 00 16 

(JeJsAHo) = (JiAM), 

where AHo =H - (H ) and AM =M - (M). Using 
this identity, we express (JeJs) in terms of the largest 
eigenvalue of the QTM A(A , Ai): 



(JeJs) = -jpd 2 Xi J d h In A(A , Ai)dA 



(6) 



A o = l,Ai=0 



To simplify the above equation further, we take the 
derivative with respect to h and Ao of both sides of eq. (4) 
after shifting the variable x — > x + i. Comparing the re- 
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sultant equations, we find 

j {d h \n V (x + i) - P)8 Xo ln(l + r]-\x + i))dx 

= <j> (d Xo lnr](x + i) - (3e(x + i))d h ln(l + ^(x + i))dx. 

Therefore we quickly see 

j> d Xa ln(l + n~ x {x + i))dx = -2irA j^lnA - || . 

Here we have used e(x) = —2irAai(x) and eq. (3). Finally 
substituting the above equation into eq. (6), we arrive at 

{JeJs) = £ ^ in(i + v ~ 1{x + i))dx - 

In Fig. 2, we show the temperature dependence of the 
current correlation (JeJs) for the isotropic point A = 1. 
As mentioned above, this correlation is strictly zero at 
h = 0. At finite magnetic fields, we find that (JeJs) has 
a finite temperature maximum which increases with in- 
creasing magnetic field for < h < h c = 4. For h > h c , 
the maximum decreases and the corresponding temper- 
ature T shifts to higher values with increasing magnetic 
field. At h — > oo in which all spins point up, the tem- 
perature To moves to infinity for fixed J. Consequently 
(JeJs) = at h = oo. For h weaker than the critical 
field h c , (JeJs) is linear in T at low-temperature. On 
the other hand, (JeJs) exponentially decays for h > h c 
because of the existence of the mass gap. 
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Fig. 2. Temperature dependence of the correlation (J7ei7s} of the 
isotropic case A = 1 for various magnetic fields. 



Low-temperature asymptotics (T <C h) .- Next we 
would like to discuss the leading contributions at low 
temperature T< 1 and T <C h (i.e. fih 3> 1) in which 
the "logarithmic correction" terms especially for the 
isotropic case are next-leading contributions. In this case 
the auxiliary function ?7 -1 (a;) is exponentially small on 
the upper contour, i.e. 1 + 2i) ~ exp(— /?M) (S > 0). 
Moreover the function on the lower contour Im x = be- 
haves as ln(l + rj^ 1 (x)) 3> 1 for the region x G [Al, Ar] 
and ln(l + rj^ 1 (x)) <C 1 for x < A L and x > Ar. There- 



fore the NLIE (4) reduces to In t](x) = (3e(x) 
where e(x) is the dressed energy given by 



0(T), 



e(x) = (\ + X 1 Ad x )e(x) + h- / a 2 (x-y)e(y)dy. (7) 

J Al 

The right (left) "Fermi point" Ar (Al) depending on 
Ao and Ai is obtained from the condition e(Ar) = 
£(Al) = 0. To consider the 0(T) correction to In rj(x), 
we take into account the behavior near the Fermi point: 
\nrj(x) = /3e'(A R/L )(x- Ar/ l ) for x ~ A R/L . Using this, 
we evaluate the 0(T) contribution 




y)Hl + V - 1 (y)) 



dy)a 2 (x- y) ln(l + ??(?/)) (8) 



by changing the variable as In r)(y) — > z and using the 
values ln?7(AR/L) = and ln?7 _1 (0) = lnr/(±oo) = oo. 
The result is 



a 2 (x-A R ) a 2 (x-A L ) 



£'(Ar) £'(A L ) 

Thus we obtain the auxiliary function In r](x) up to 0(T): 

\nn(x) - 3e(x) + — f a ^ x ~ A ^ _ a ^ x Z A ^ ' 
V[) ma;j+ 6/3\ s'(Ar) e '(A L ) 

Applying the same procedure and using the above ob- 
tained relations, we have the low-temperature behavior 
of the eigenvalue (3) up to 0(T): 

f-A R 

a\{x)e{x)dx 



InA=^-/3 
2 



A L 



+ 6/3\ e '(AR) e '(A L )J + U(1 J ' (9J 
where p(x) is the density function given by 

/•Ar. 

p(x) = a 1 (x)- a 2 (x - y)p{y)dy. (10) 

JAl 

Combining (9) with (5) and (6), we obtain and 
(JeJs) for f3h > 1: 

(ttA) 2 T 3 / a'(A)p(A) , a(A)^(A) 

e' 2 (A) 



(Je 2 ) = ^s(0)T + 



2a(A)e"(A)/o(A) 



e' 2 (A) 
a(A) 



+ 



J (A)p(A)\ 



(JeJs) =hD s (0)T + 



e ,3 (A) 27rAe'(A) e' 3 (A) J ' 

7rAr 3 [ a 2 (A)£(A) | a'(A)g(A) 



6 1 e' 3 (A) 



e' 2 (A) 
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+ 



<*(A)g'(A) 
(A) 



r'2 



2a(A) £ "(A)g(A) ) 
e' 3 (A) J' 



(11) 



where A is the Fermi point determined by e(±A) = 
(note that we set (A , Ai) = (1, 0) in (7) and (10)); D s (0) 
is the zero temperature spin stiffness D s (0) — vf£, 2 (A)/it; 
vf is the Fermi velocity defined by vf = e'(A)/(27rp(A)); 
the function a(x) and the dressed charge £(x) are given 

by 



a(x) — e"(v) 



a 2 (x - y)a(y)dy, 



a2(x - y)£{y)dy. 



Drude weights and Seebeck Coefficient.- Finally we dis- 
cuss the thermal Drude weight D t h (T) and the magnetic 
Seebeck coefficient S at low-temperature {j3h 3> 1). To 
evaluate them, we have to explicitly determine the spin 
Drude weight D S (T) as well as the above obtained cur- 
rent correlations and (JfJs)- Though the Drude 
weight for the spin transport at finite temperatures and 
zero magnetic field has been evaluated by Zotos already 
a while ago, 7 the validity of the results is still debated. 
Here, we determine the low-temperature behavior alter- 
natively without directly calculating D S (T). Utilizing the 
phcnomcnological relation between the thermal conduc- 
tivity k and the specific heat C, i.e. k = C(T)v f t (t: 
relaxation time) , we assume the low-temperature asymp- 
totics of the thermal Drude weight as 



D th (T) = ^T + 0(T 2 ). 



(12) 



Here we have used k = D th (T)r and C = 7rT/(3v F ). 20 
Note that t = oo for the present case. The validity of 
eq.(12) is verified at h = 11 and for the XY model (A = 
0) for < h < h c (h c = 2). Combining (12) with (1) and 
using the result (11), we determine the low-temperature 
behavior of D B (T). From the resultant equation and (11), 
we obtain the ratio of {JfJs) and D S (T); 

ttT 3 



(JeJs) 



D S (T) QAp(A)aA) \ 



+ 



Aa(A) 
27rp(A)uj 



Finally substituting this into (2), we arrive at the lead- 
ing low-temperature behavior of the magnetic Seebeck 
coefficient for T <C h: 



S(T) 



ttT 



1 



Aa(A) 



0(T 2 ). 



6Ap(A)£(A) { 2np(A)v F 

In Fig. 3, the coefficient of the leading correction is de- 
picted as a function of the magnetic field for various 
anisotropy parameters. For weak interaction strengths 
A < A ~ 0.5, the leading behavior is negative. On 
the contrary, for A > A , due to effects of strong cor- 
relations, the Seebeck coefficient changes sign at certain 
magnetic field ho (or equivalcntly certain magnetization 
Mo)- The value ho shifts to higher values with increasing 
the interaction strengths. Recently similar behavior has 
also been observed in the electric Seebeck coefficient of 



the Hubbard model for finite U and T by the numerical 
diagonalization for small systems. 21 

In summary, we have discussed the magnetothermal ef- 

1.5 



CO 




Fig. 3. The leading contribution to the Seebeck coefficient at low 
temperatures T -C h. 



fects and the thermal transport at finite magnetic fields. 
The magnetic Seebeck coefficient changes sign above cer- 
tain interaction strengths, because of effects of strong 
correlations. 
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